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We show that the numerical renormalization group (NRG) is a viable multi-band impurity solver 
for Dynamical Mean Field Theory (DMFT), offering unprecedent real-frequency spectral resolution 
at arbitrarily low energies and temperatures. We use it to obtain a numerically exact DMFT solution 
to the Flund’s metal problem for a three-orbital model with filling factor na = 2. The ground state 
is a Fermi liquid. The one-particle spectral function shows a coherence-incoherence crossover with 
increasing temperature, with spectral weight being transfered from low to high energies, and evolves 
qualitatively differently from a doped Mott insulator. The spectral function has a strong particle-hole 
asymmetry, and in the incoherent regime the one-particle self-energy shows approximate power- 
law behavior for positive frequencies only. The spin and orbital spectral functions show “spin- 
orbital separation”: spin screening occurs at much lower energies than orbital screening. The 
renormalization group flows clearly reveal the relevant physics at all energy scales. 


PACS numbers: 71.27.+a, 71.10.Fd, 75.20.Hr 

Introduction .— A widely-used method for dealing 
with interactions in strongly-correlated electron systems 
and electronic structure calculations is dynamical mean 
field theory (DMFT) [1, 2]. It treats the interplay be¬ 
tween a given lattice site (the “impurity”) and the rest 
of the lattice (the “bath”) as a quantum impurity model 
with a self-consistently determined hybridization func¬ 
tion. Since DMFT’s performance depends on that of the 
method used to solve this impurity model, much effort 
has been invested over the years to develop ever more 
powerful impurity solvers. For multi-band models, con¬ 
tinuous-time Quantum Monte Carlo (ctQMC) methods 
appear to be the current favorites in terms of versatility 
and performance [3]. However, they are not without lim¬ 
itations: sign problems can occur, low-temperature cal¬ 
culations are costly, and obtaining real-frequency spec¬ 
tra requires analytic continuation of imaginary (Matsub- 
ara) frequency QMC data, which is notoriously difficult. 
Thus, there is a continued need for real-frequency impur¬ 
ity solvers suitable for multi-band DMFT applications. 

In this work, we show that the numerical renormaliza¬ 
tion group (NRG) [4-6] is such a tool, offering unpreced¬ 
ented real-frequency spectral resolution at low energies. 
NRG is the gold standard for impurity models, with nu¬ 
merous previous DMFT applications (e.g. [7-13]), but 
until recently was limited to models with at most two 
bands. However, recent technical progress [14-16] has 
now made three-band calculations feasible [17-19]. 

We illustrate the potential of DMFT-I-NRG by study¬ 
ing a three-band model [20-22] with both a Hubbard in¬ 
teraction U and a ferromagnetic Hund’s coupling J, with 
U(l)^j^ X SU(2),,p X SU(3)^,.^ symmetry for its charge (ch), 
spin (sp) and orbital (orb) degrees of freedom. This is 
the simplest model of a 3-band “Hund’s metal” [23, 24]. 
These are multi-orbital materials with broad bands which 
are correlated via the Hund-J rather than the Hubbard-U 
interaction. Examples are iron pnictide and chalcogen- 


ide high-temperature superconductors [23, 25], ruthen- 
ates [26, 27], and other 4d transition metal oxides [21, 28]. 

Early DMET studies using continuous-time QMC 
(ctQMC) [3] as impurity solver suggest that consequences 
of the Hund’s rule coupling include (i) Fermi-liquid beha¬ 
vior at low energies [23] and (ii) a coherence-incoherence 
crossover with increasing temperature [23], relevant for 
various material systems [27, 29]. The incoherent regime 
is characterized by (iii) fractional power laws [20, 26, 30] 
and (iv) the coexistence of fast quantum mechanical or¬ 
bital fluctuations and slow spin fluctuations [20]. How¬ 
ever, since ctQMC can not reach truly low temperatures, 
(i) could not be conclusively established yet, and a more 
detailed understanding of (ii-iv) is difficult to achieve 
based on imaginary-frequency data alone. Our real- 
frequency DMFT-I-NRG results definitively settle these 
issues and yield further insights. We find (i) a Fermi 
liquid ground state; a real-frequency one-particle spec¬ 
tral function showing (ii) a coherence-decoherence cros¬ 
sover with significant transfer of spectral weight from 
low to high energies, and (iii) strong particle-hole asym¬ 
metry, which leads to the above-mentioned apparent frac¬ 
tional power laws; and (iv) two-stage screening, where 
spin screening occurs at much lower energies than orbital 
screening (“spin-orbital separation”), (v) The T = 0 
spectral properties are similar with or without DMFT 
self-consistency, in marked contrast to Mott-Hubbard 
systems, where the DMFT self-consistency opens a gap in 
the quasiparticle spectrum at large interaction strength. 

Model :— Our three-band model has Hamiltonian 

^ = + X (la) 

i {ij)u 

Hintldl] = I JiV, + ^{u- \J) W(W - 1) - J^l (lb) 

Here creates an electron on site i of flavor (fl) v, with 
composite index v = (ma) labelling its spin {a =t,4.) 
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and orbital (m = 1,2,3). Ni = is the total 

number operator for site i and its total spin, with 
components = J^maa' d\,na\'^aa'dima', where (t“ are 
Pauli matrices. We study a Bethe lattice with nearest- 
neighbor hopping amplitude t, used as energy unit (t = 
1). Onsite interactions are described by l?int. The onsite 
Coulomb interaction U penalizes double occupancy. The 
ferromagnetic coupling J > 0 accounts for Hund’s first 
rule by favoring a large spin per site. We choose the 
chemical potential /i such that the filling per lattice site 
is one below half-Hlling, (Ni) ~ 2, inducive to an intricate 
interplay of spin and orbital degrees of freedom. 

Methods. — We use single-site DMFT to map the lat¬ 
tice model onto a 3-band Anderson-Hund model (AHM) 
of the form Harm = SdN + Hint[dl] + ffbath- Here dl 
creates a local (“impurity”) electron of flavor i/ with en¬ 
ergy £d = —/i, experiencing local interactions H-mt, with 
total number and spin operators N and S defined analog¬ 
ously to Ni and S^. The local site on average hosts two 
electrons (n^ = (N) ~ 2), forming a spin triplet and or¬ 
bital triplet (the one hole relative to half-hlling can be in 
one of three orbital levels). The local electrons hybridize 
with a 3-band spinful bath, 

.Hbath = ’ (2) 

ku 

with a hybridization function r(e) = tt \Vk\^6{e — Sk) 
that fully characterizes the impurity-bath interplay. In 
DMFT, F(e) has the role of the effective Weiss mean 
field and is determined self-consistently [1, 2, 31]. We 
studied both the self-consistent AHM (scAHM), and for 
comparison also the pure impurity AHM (iAHM) without 
self-consistency, using a flat density of states with half¬ 
bandwidth D, r(£) = r0(D — jej). 

We use full-density-matrix (fdm) NRG [32] exploiting 
non-Abelian symmetries [15], both to solve the iAHM 
and for each scAHM iteration (for NRG details, see [31]). 
The key idea of NRG, due to Wilson [4], is to discret¬ 
ize the bath’s continuous spectrum logarithmically, map 
the model onto a semi-infinite “Wilson” chain with ex¬ 
ponentially decaying hopping amplitudes, and exploit 
this energy-scale separation to iteratively diagonalize the 
model while discarding high-energy states. This allows 
one to zoom in on low-energy properties, at the expense 
of having only coarse-grained resolution at high energies. 
Nevertheless, NRG results are accurate also for spectral 
integrals even if these include large energies, since they 
can be evaluated using discrete, unbroadened NRG data. 

Matsubara benchmark. — We illustrate this by bench¬ 
marking NRG versus ctQMC [31], which treats the bath 
as a continuum and has no bath discretization issues. We 
used both methods to compute the self-energy S(iwn) 
of the Matsubara correlator G(iaj„) associated with the 
retarded local correlator G^{oj) = {d^\\dl)uj. In NRG, 
its spectral function is expressed in terms of discrete 


scAHM : U = 5,J = 1, fi = 7.032, T = 0.001 
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Figure 1. (Color online) Benchmark comparison of NRG and 
ctQMC for the 3-band scAHM. (a) Imaginary and (b) real 
part of the self-consistently converged self-energy as function 
of Matsubara frequencies. Grey lines in (a) are power-law fits 
to low and intermediate frequency data, respectively. The 
inset of (a) shows ImE(ica„) on a linear scale. 


data, A{uj) = Ss®s<5(w — Cs), hence 

G(ia;„) = / da;A(a;)/(iw„ - w) - ^s)- 

Fig. 1 compares NRG and ctQMC results for E(ia;„) at 
T = 0.001. The agreement is excellent, also at large 
frequencies. The numerical costs differ vastly, however, 
~ 10^ vs. 10® CPU hours [31], since the chosen temper¬ 
ature is challengingly low for ctQMC, whereas NRG can 
access any temperature. Luttinger pinning at zero fre¬ 
quency [1, 33] is fulfilled within 1.5% for both methods 
[31]. ImE(icij„) displays fractional power-law behavior 
for intermediate frequencies (0.1 ^ < 1), as found in 

[20, 26], and Fermi-liquid behavior for very low frequen¬ 
cies, ImE(ia;„) oc uin, as found in [23, 27]. 

Coherence-incoherence crossover. We now turn to 
real-frequency properties [Fig. 2]. At zero temperat¬ 
ure, the local spectral function A{uj) of the scAHM 
shows a well-defined low-energy quasiparticle peak and 
—ImE^(a;) a dip reaching down to zero [insets of 
Figs. 2(a,b)]. This indicates that strong Kondo-type 
screening correlations exist between bath and local spin 
and orbital degrees of freedom. At higher energies, A{ijj) 
also shows incoherent, rather flat particle-hole asymmet¬ 
ric side peaks, that reflect charge fluctuations. With 
increasing temperature, a coherence-incoherence cros¬ 
sover occurs: the quasiparticle peak weakens and even¬ 
tually gives way to a pseudogap [Fig. 2(a)], and con¬ 
currently the dip in —lmS^(a;) is smeared out into a 
broader minimum, which eventually evolves into a max¬ 
imum [Fig. 2(b)]. During this crossover, quasiparticle 
weight is transferred from low to high energies, in a 
way reminiscent of recent photoemission measurements 
[34] (though the interpretation offered there invoked an 
orbital-selective Mott phase). 

The behavior seen in Figs. 2(a,b) is characteristic of 
“Hund’s physics”, arising when the Hund’s coupling J is 
sizeable {J/U nonzero and not too small). It contrasts 
markedly with the behavior associated with a Mott trans¬ 
ition, as exemplified, e.g., in Figs. 2 and 3 of [9] for the 
1-band Hubbard model: the side peaks in our A[(jj) are 
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Figure 2. (Color online) (a) The local spectral function 
A{iu) and (b) the imaginary part of the retarded self-energy, 
ImE^(a;), for the scAHM, plotted versus frequency for four 
temperatures. Insets show a larger frequency range for T = 
10“®. (c,d) Same as in (a,b), but for an iAHM. 

rather flat, the maxima in our —ImE^(w) never diverge, 
and the spectral weight near a; ~ 0 remains nonzero at all 
temperatures, implying that metallic behavior persists. 
In this sense, Hund’s metals exhibit properties quite dis¬ 
tinct from Mott-Hubbard systems. 

Since the scAHM is based on an impurity model, it is 
instructive to study a corresponding iAHM, with para¬ 
meters tuned to yield a similar spectral function at T = 0 
[Fig. 2(c,d)]. It likewise features a large low-energy 
(Kondo) peak that weakens with increasing temperat¬ 
ure, though no pseudogap occurs. This shows that the 
T = 0 spectral properties of the scAHM are governed 
by the impurity physics of the iAHM, while the scAHM 
transfer of spectral weight with increasing T is brought 
about by DMFT self-consistency. 

Particle-hole asymmetry. We next exploit the power 
of NRG to zoom in to arbitrarily low energy scales: In 
Figs. 3(a,b), we replot on a logarithmic scale the data 
(black/red for sc/iAHM) from Figs. 2(a-d) for Aiuj) and 
Im at T = 10“®. For comparison, the right column 

of Fig. 3 again shows results for the iAHM, but using 
parameters that yield smaller crossover scales (defined 
below), to better separate the low-energy features asso¬ 
ciated with spin and orbital screening from high-energy 
features associated with charge fluctuations. Note again 
the striking qualitative similarity between the scAHM 
(black) and iAHM (red/blue) spectra - clearly, for T ~ 0 
DMFT self-consistency plays no major role. 

With decreasing temperature, the quasiparticle peaks 
in Figs. 2(a,c) show an increasing particle-hole asym¬ 
metry, not surprising away from half-filling, which at 
T ~ 0 is very pronounced: in Figs. 3(a,d) for A(a;), the 
thick (w < 0) lines show a shoulder-like structure for 



Figure 3. (Color online) (a-f) Spin-orbital separation in 
real-frequency, ground state correlators. The left column 
uses the same parameters and color code as Fig. 2 for the 
scAHM (black) and iAHM with (r,D) = (0.778,0.5) (red). 
For comparison, the right column shows iAHM results with 
(r, D) = (0.200,1.0) (blue), yielding smaller crossover scales. 
(a,d) The local spectral functions, (b,e) the local self-energy, 
and (c,f) the spin and orbital susceptibilities, Xsp (solid) and 
Xorb (dashed). We use a logarithmic frequency scale, with 
thick (thin) lines for tj < 0 (a; > 0). Insets show data on a 
linear scale. In all panels, solid (dashed) vertical lines mark 
the spin (orbital) Kondo scale, Grey guide-to-the 

eye lines indicate Fermi-liquid power laws (solid) or apparent 
fractional power laws (dashed). Inset to (f): Kondo scales 
(solid) and (dashed) for the iAHM, plotted as function 
of Ud- (g,h) NRG eigenlevel flow diagrams for the scAHM and 
iAHM of panels (a-c) and (b-d), respectively, showing the res¬ 
caled energies of the lowest-lying eigenmultiplets of a Wilson 
chain of (even) length k, plotted versus its characteristic level 
spacing Uk = (see text). Numbers above lines give 

multiplet degeneracies, Q their symmetry labels. 

intermediate frequencies (between the vertical solid and 
dashed lines), while the thin (w > 0) lines do not; and in 
Figs. 3(b,e) for ImE'^(a;), the thick lines show a plateau¬ 
like structure, whereas the thin lines show approximate 
~ power-law behavior (with non-universal a). For 
the Matsubara self-energy obtained via the Hilbert trans¬ 
form E(itLin) = —^ f dujim (lu) / (iujn — w), the asym¬ 
metric contributions from the power-law and shoulder 
in ImE^(a; < 0) conspire in such a way that ImE(iaj„) 
shows an apparent fractional power law, but ReE(iaj„) 
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does not [Fig. l(a,b)]. Conversely, this example ilustrates 
that care is due when drawing real-frequency conclusions 
from imaginary-frequency power laws: if one is present 
only for Iml](iw„), not ReS(iw„) (as in [20, 26]), then 
Iml]^(w) need not show pure power-law behavior. 

Spin-orbital separation. Next we elucidate the screen¬ 
ing of local spin and orbital degrees of freedom by the 
bath of conduction electrons. To this end. Figs. 3(c,f) 
respectively show the imaginary part (%") of the dynam¬ 
ical spin and orbital susceptibility, Xsp = 
and Xorb = orbital operators T“ = 

Emm'crdLA'^^m'dm’a, where are the SU(3) Gell- 
Mann matrices, normalized as Tr[r“r^] = 2Sab- Both x"p 
and x ",.|3 exhibit a peak with (nearly) power-law flanks, 
characteristic of Kondo screening of the local spin and 
orbital degrees of freedom. Strikingly, for both scAHM 
and iAHM the peak for x"p occurs at a much lower en¬ 
ergy and is much higher than for Xorb- peak posi¬ 
tions lie near the spin and orbital Kondo scales, and 
(vertical solid and dashed lines), calculated from 
Xsp(w) and Xorb(<j-’) using the recipe of [18]. acts as 
the coherence scale below which Fermi-liquid behavior 
(ImS^(a;) oc x"p,orb see Figs. 3(b-f), grey lines) 
sets in. The SU(2)sp and SU(3)orb crossover scales dif¬ 
fer strongly, , because the Kondo temperature 

for an SU(A^) Kondo model scales as InTx — 1/iV. This 

implies two-stage screening, with spin screening occur- 
ing at significantly lower energies than orbital screening. 
This “spin-orbital separation”, featuring a very small co¬ 
herence scale and an intermediate regime with screened 
orbital degrees of freedom coupled to slowly fluctuating, 
large spins, is a central result of our work. 

The inset of Fig. 3(f) depicts and for the 
iAHM as function of the filling For ~ 1, where 
the bare impurity’s ground state has SU( 6 ) symmetry 
also for J ^ 0, ~ ~ As Ud increases 

from 1 to 2 , and split apart and spin-orbital 
separation sets in. As continues to increase towards 
3, becomes very large (> D), reflecting the fact that 
for half-hlling the orbitals form an orbital singlet from the 
outset. In this sense, ~ 2 is special: there conditions 
are optimal for the Hund coupling to align two spins in 
different orbitals without forming an orbital singlet. 

RG flow. In renormalization group (RG) terms, the 
two-stage screening discussed above is associated with 
the RG flow between three fixed points, describing high-, 
intermediate- and low-energy energy excitations. Their 
effective fixed point Hamiltonians have ground state 
multiplets whose spinxorbit structure is tripletxtriplet, 
tripletXsinglet and singletxsinglet, implying a ground 
state entropy of ln(9), ln(3) and ln(l), respectively. The 
RG flow between these fixed points can be visualized via 
NRG eigenlevel flow diagrams [Fig. 3(g,h)]. Technically, 
they show how the lowest-lying rescaled eigenlevels of a 
length-fc Wilson chain evolve with k, where “rescaled” 


means given in units of Wfc = (A > 1 being a dis¬ 

cretization parameter). Conceptually, these levels rep¬ 
resent the finite-size spectrum of the impurity-|-bath put 
in a spherical box of radius ~ A^/^, centered on the 
impurity [4, 35]: as k increases, the finite-size level spa¬ 
cing LVk (x. l/Rk decreases exponentially. The structure of 
the finite-size spectrum is fixed (/c-independent) while Wfc 
lies within an energy regime governed by one of the fixed 
points, but changes when passes through a crossover 
scale between two fixed points. 

Figs. 3(g,h) show this RG flow for the scAHM and the 
iAHM, revealing similar behavior for both [36]. We label 
multiplets by their U(l)^jj x SU(2)gp x SU(3 )q,.j, symmetry 
labels, Q = [g, 25', < 71 ^ 2 ]; here q denotes particle number 
relative to half-filling, S spin, and (( 7192 ) an SU(3) irrep, 
identified by a Young diagram with qi (72 (( 72 ) boxes in 
its first (second) row. The flow of the lowest-lying levels 
reveals two crossover scales, and (whose spa¬ 

cing, though, is too small for the level flow in between to 
become stationary, i.e. fc-independent). As uik drops be¬ 
low T^'^, orbital screening sets in, favoring orbital sing¬ 
lets [( 51 ( 72 ) = (00)], hence other multiplets rise in energy, 
the more so the larger ( 5152 )- Similarly, as uJk drops be¬ 
low spin screening sets in, favoring spin singlets and 
pushing up multiplets with 5 7 ^ 0. For <C T^, the 
ground state is a spin and orbital singlet [Q = (0, 0, 00)]. 
We have checked that its excitation spectrum can be in¬ 
terpreted in terms of non-interacting single-particle ex¬ 
citations, thus confirming its Fermi-liquid nature. 

Conclusions. We have demonstrated the potential of 
DMFT-I-NRG as real-frequency method to treat multi¬ 
orbital systems, with no need for analytic continuation. 
Applied to a Hund’s metal, it revealed subtle spectral fea¬ 
tures which are manifestly different from those of Mott- 
Hubbard systems, and which can be probed in photoe¬ 
mission and STM spectroscopies. Our work is a first step 
towards using LDA-I-DMFT-I-NRG to calculate AG and 
DC transport properties in strongly correlated materi¬ 
als. Such applications will typically involve less orbital 
symmetries than the model studied here, but could be 
treated using the recent “interleaved” NRG approach of 
[16]. The latter yields results of comparable accuracy and 
efficiency as when symmetries can be exploited [37] , and 
calculations with up to 5 bands seem within reach. 

A key advantage of NRG is its ability to iteratively un¬ 
cover the system’s RG flow from high to low energies, re¬ 
vealing the relevant physics at each scale. In the context 
of Mott-Hubbard systems, RG ideas have been very fruit¬ 
ful even in very approximate implementations [38-40]. 
For the present Hund’s metal, the numerically exact RG 
flow achieved via DMFT-I-NRG revealed a clear, simple 
picture of the crossover from the incoherent to the coher¬ 
ent Fermi-liquid regime: two-stage screening of first or¬ 
bital, then spin degrees of freedom. Using DMFT-I-NRG 
to gain this type of RG understanding of real material 
properties would be a worthwhile goal for future research. 
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SUPPLEMENTARY MATERIAL 

Recent NRG progress. While Wilson’s NRG [4, 5] has 
been tremendously successful in the past, significant fur¬ 
ther progress was achieved in recent years, triggered by 
the realizations that it (i) can be formulated in MPS lan¬ 
guage [6, 42], and (ii) that the discarded states can be 
used to construct a complete many-body basis of approx¬ 
imate energy eigenstates [43] . It has become possible (iii) 
to reliably calculate finite-temperature spectral functions 
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in sum-rule conserving fashion [32, 44]; (iv) to treat a 
bath with nontrivial hybridization function by suitably 
optimizing its representation in terms of discrete bath 
states [45-47]; to treat multi-band models, either (v) by 
exploiting non-Abelian symmetries [14, 15] if the bands 
couple symmetrically to the impurity [17-19], as here, 
or (vi) by using an “interleaved” discretization scheme 
[16, 37]; and (vii) to greatly increase numerical efficiency 
by optimizing MPS bond dimensions and to ensure ac¬ 
curacy by checking the discarded weight [48]. Taken to¬ 
gether, these advances make NRG a highly competitive 
real-frequency impurity solver for DMFT. 

Details of our NRG implementation. Wilson’s NRG 
approach is based on logarithmically coarse-graining the 
continuous bath in energy space into intervals of expo¬ 
nentially decreasing widths in order to resolve even the 
lowest relevant energy scale of the impurity system. To 
accurately represent the frequency-dependent hybridiza¬ 
tion function of the DMFT self-consistency condition in 
terms of a set of discrete states, we use a numerically 
stable implementation [47] of the discretization scheme 
of Ref. [45, 46]. 

The model is then mapped onto a semi-infinite 
“Wilson” chain with exponentially decaying hopping 
amplitudes. This energy-scale separation is exploited 
to iteratively diagonalize the model by adding one site 
at a time, while discarding high-energy states. The ac¬ 
curacy of this truncation procedure can be checked by 
estimating the discarded weight (5pdisc [48], a quantitat¬ 
ive convergence criterion inspired by DMRG, after each 
run. Empirically, an NRG run is well-converged when 
Spdisc < 10 “^^. 

We exploit all available non-Ablian symmetries of the 
AHM studied here, namely U(l),,.j, x SU(2)gp x SU(3 )q,.j,, 
by using the QSpace approach developed by A. Weichsel- 
baum [15]. QSpace is a tensor library that is able to treat 
Abelian and non-Abelian symmetries on a generic level. 
In the presence of symmetries, the state space can be 
organized into symmetry multiplets, and tensors “factor¬ 
ize” into two parts, acting in the reduced multiplet space 
and the Glebsch Gordon coefficient space, respectively, 
vastly reducing numerical cost. 

We use full-density-matrix (fdm)NRG [6, 32] to calcu¬ 
late high-quality, sum-rule-conserving spectral functions 
at arbitrary temperature. To implement the DMFT self- 
consistency loop (described below) we smoothen the dis¬ 
crete spectral data provided by fdmNRG using the log- 
Gaussian smoothening approach of Ref. [32] . To improve 
the resolution of the smoothened spectral data we calcu¬ 
late the self-energy in every iteration from the ratio of two 
correlation functions [8]. Moreover, for the last DMFT 
iteration we average over two versions of the discretiz¬ 
ation grid (^-averaging with W = 2, [45, 46, 49]). We 
have explored also doing z-averaging for earlier DMFT- 


iterations, but that did not noticably improve conver¬ 
gence. 

NRG-related computational parameters. The per¬ 
formance of NRG is governed by the following compu¬ 
tational parameters: the dimensionless logarithmic dis¬ 
cretization parameter A; the truncation energy Atrunc in 
rescaled units (as defined in Ref. 6), up to which all eigen- 
multiplets are kept, unless their number exceeds 
the maximal number of kept multiplets per iteration; 
the number of z-shifts for z-averaging ; and the log- 
Gaussian broadening parameter a for smooth spectral 
data [32] . 

We use A = 4, cr = 0.8 and Atrunc = 7. For impur¬ 
ity spectral functions and self-energies we used = 

2500, yielding discarded weights of Spdisc < 10“^^. For 
the susceptibilities Xsp and Xorb we use smaller value for 
iV^ggp = 1500 (because the operators involved are numer¬ 
ically rather costly); this yields values for and 
that are consistent with the crossover scales derived from 
NRG eigenlevel flow diagrams. 

DMFT self-consistency loop. — In single-site DMFT 
a quantum lattice model is treated in a quantum mean- 
field fashion. Spatial correlations are frozen out by 
dropping the momentum-dependence of the lattice self¬ 
energy, whereas temporal quantum fluctuations are re¬ 
tained. The lattice dynamics is then fully captured 
by the local retarded lattice Green’s function, Gll^^{uj), 
given in terms of the purely local but still frequency de¬ 
pendent self-energy, - or equivalently by the re¬ 

tarded Green’s function of an impurity model G(^p{uj) 
with equal local interactions (equal £-^( 0 ;)) and effect¬ 
ive hybridization r(a;). This equivalence, Gj^^^(a;) = 
Gimp(w) = G(cli), imposes a self-concistency condition, 
that simplihes to r(a;) = —t^ImG;(^p(w) for a Bethe lat¬ 
tice, and fully maps the quantum lattice problem onto 
an effective quantum impurity problem by iteratively 
determining r(a;) [1]. In each iteration of the DMFT 
self-consistent mapping, we solve the quantum impurity 
model with NRG. 

Gomputational costs. Our ctQMC solver, developed 
by K. Haule, is based on an expansion in the hybridiza¬ 
tion function [50, 51]. Its implementation is described in 
[52]. For the benchmark calculations shown in Fig. I, our 
ctQMC solver was run in parallel using 2500 cores with 
10® Monte Carlo steps in each core. To achieve DMFT 
self-consistency for T = 0.001, it needed 20 iterations of 
6 hours each. Our NRG code, developed by A. Weichsel- 
baum [6, 15], achieved DMFT self-consistency after 7 it¬ 
erations of 10 hours each, on a single 8-core machine with 
128 GB memory. This results in nearly three orders of 
magnitude better numerical efficiency. Moreover, NRG 
is able to access arbitrarily low temperatures, whereas 
for ctQMC the numerical costs grow exponentially with 
decreasing temperature. 


